Local Ensemble Transform Kalman Filter Research Articles

Impact of atmospheric forcing on SST biases in the LETKF-based ocean research analysis (LORA)

In the previous study, the authors have produced an eddy-resolving ocean ensemble analysis product called the local ensemble transform Kalman filter (LETKF)-based ocean research analysis (LORA) over the western North Pacific and Maritime Continent regions using an ocean data assimilation system driven by the Japanese operational atmospheric reanalysis dataset known as the JRA-55. However, the LORA includes warm biases in sea surface temperatures (SSTs) in coastal regions during the boreal winter. In this study, we perform sensitivity experiments with atmospheric forcing using an ocean forcing dataset known as the JRA55-do, which adjusts the JRA-55 to high-quality reference datasets to reduce biases and uncertainties. The results show that the nearshore warm SST biases are significantly improved by the JRA55-do. During the boreal autumn, the improvement comes from mainly two factors: (i) enhancement of surface cooling by latent heat releases caused by removing contamination of weak winds at the land grid cells, and (ii) weakening surface heating by downward shortwave radiation through the adjustment in the JRA55-do.During the boreal winter, enhanced cooling by the analysis increments suppresses the growth of the warm SST biases when the JRA55-do is used. However, if the JRA-55 dataset is used, the adaptive observation error inflation (AOEI) scheme acts negatively to keep the nearshore SST biases in winter. Based on the innovation statistics, the AOEI inflates the observation errors when the differences between the squared observation-minus-forecast innovations and the squared forecast ensemble spreads are larger than the prescribed observation error variance, and improves the accuracy in the open ocean, especially around the frontal regions. However, when substantial warm SST biases are formed in the previous season, AOEI's observation error inflation makes the analysis increments smaller and cannot suppress the warm biases.We also validate the analysis accuracy using various data such as sea surface height and horizontal velocities and find that the JRA55-do has significant advantages. Therefore, continuous maintenance and development of ocean forcing datasets are essential for ocean modeling and data assimilation.

Improving Afternoon Thunderstorm Prediction over Complex Terrain with the Assimilation of Dense Ground-Based Observations: Four Cases in the Taipei Basin

Abstract In this study, we investigate the impact of assimilating densely distributed Global Navigation Satellite System (GNSS) zenith total delay (ZTD) and surface station (SFC) data on the prediction of very short-term heavy rainfall associated with afternoon thunderstorm (AT) events in the Taipei basin. Under weak synoptic-scale conditions, four cases characterized by different rainfall features are chosen for investigation. Experiments are conducted with a 3-h assimilation period, followed by 3-h forecasts. Also, various experiments are performed to explore the sensitivity of AT initialization. Data assimilation experiments are conducted with a convective-scale Weather Research and Forecasting–local ensemble transform Kalman filter (WRF-LETKF) system. The results show that ZTD assimilation can provide effective moisture corrections. Assimilating SFC wind and temperature data could additionally improve the near-surface convergence and cold bias, further increasing the impact of ZTD assimilation. Frequently assimilating SFC data every 10 min provides the best forecast performance especially for rainfall intensity predictions. Such a benefit could still be identified in the earlier forecast initialized 2 h before the start of the event. Detailed analysis of a case on 22 July 2019 reveals that frequent assimilation provides initial conditions that can lead to fast vertical expansion of the convection and trigger an intense AT. This study proposes a new metric using the fraction skill score to construct an informative diagram to evaluate the location and intensity of heavy rainfall forecast and display a clear characteristic of different cases. Issues of how assimilation strategies affect the impact of ground-based observations in a convective ensemble data assimilation system and AT development are also discussed. Significance Statement In this study, we investigate the impact of frequently assimilating densely distributed ground-based observations on predicting four afternoon thunderstorm events in the Taipei basin. While assimilating GNSS-ZTD data can improve the moisture fields for initializing convection, assimilating surface station data improves the prediction of rainfall location and intensity, particularly when surface data are assimilated at a very high frequency of 10 min.

Strongly Versus Weakly Coupled Data Assimilation in Coupled Systems With Various Inter‐Compartment Interactions

AbstractCoupled data assimilation (CDA) has been attracting researchers' interests to improve Earth system modeling. The CDA methods are classified into two: weakly coupled data assimilation (wCDA), which considers cross‐compartment interaction only in a forecast phase, and strongly coupled data assimilation (sCDA), which additionally uses other compartment's information in an analysis phase. Although sCDA can theoretically provide better estimates than wCDA since sCDA fully uses inter‐compartment covariances, the potential of sCDA in practice is still in debate. We investigate conditions under which sCDA is effective by applying Local Ensemble Transform Kalman Filter (LETKF) to a coupled Lorenz96 model. Especially, we aim to generalize CDA's behavior on coupling strength basis. By continuously changing a cross‐compartment interaction in the coupled Lorenz96 model, we find that sCDA's superiority over wCDA is particularly evident with large cross‐compartment interactions which attenuates dynamics of each compartment (the increase of variables in one compartment induces the reduction of those in the other compartment) and does not intensify the chaoticity of the whole system. The performance of sCDA is more sensitive to the LETKF's hyperparameters (i.e., localization and inflation) than wCDA. Furthermore, through two imperfect model experiments in which model parameters are uncertain, we quantify sCDA's vulnerability to model inaccuracies. The findings of the relationship between the skill of sCDA and inter‐compartment interactions give implications for future CDA strategies.

Stochastic perturbations in the semi-Lagrangian advection algorithm of the SL-AV global atmosphere model

Abstract An algorithm for stochastic perturbation of the semi-Lagrangian trajectories is implemented in the ensemble weather prediction system based on the global atmosphere model SL-AV20 with a horizontal resolution of approximately 20 km, 51 vertical levels, and Local Ensemble Transform Kalman Filter (LETKF). The combined use of methods for stochastic perturbation of trajectories and the parameters and tendencies of subgrid-scale processes parameterizations allows to generate ensembles with a larger spread compared to ensembles without stochastic perturbations of trajectories. An improvement in probabilistic estimates of the ensemble forecasts for various variables is shown. The comparison of two versions of ensemble prediction system is presented.

Local Volume Solvers for Earth System Data Assimilation: Implementation in the Framework for Joint Effort for Data Assimilation Integration

AbstractThe Joint Effort for Data assimilation Integration (JEDI) is an international collaboration aimed at developing an open software ecosystem for model agnostic data assimilation. This paper considers implementation of the model‐agnostic family of the local volume solvers in the JEDI framework. The implemented solvers include the Local Ensemble Transform Kalman Filter (LETKF), the Gain form of the Ensemble Transform Kalman Filter (GETKF), and the optimal interpolation variant of the LETKF (LETKF‐OI). This paper documents the implementation strategy for the family of the local volume solvers within the JEDI framework. We also document an expansive set of localization approaches that includes generic distance‐based localization, localization based on modulated ensemble products, and localizations specific to ocean (based on the Rossby radius of deformation), and land (based on the terrain difference between observation and model grid point). Finally, we apply the developed solvers in a limited set of experiments, including single‐observation experiments in atmosphere and ocean, and cycling experiments for the atmosphere, ocean, land, and aerosol assimilation. We also illustrate how JEDI Ensemble Kalman Filter solvers can be used in a strongly coupled framework using the interface solver approximation, which provides increments to the ocean based on observations from the ocean and atmosphere.

Assessing the assimilation of Himawari-8 observations on aerosol forecasts and radiative effects during pollution transport from South Asia to the Tibetan Plateau

Abstract. Emissions from South Asia (SA) represent a critical source of aerosols on the Tibetan Plateau (TP), and aerosols can significantly reduce the surface solar energy. To enhance the precision of aerosol forecasting and its radiative effects in SA and the TP, we employed a four-dimensional local ensemble transform Kalman filter (4D-LETKF) aerosol data assimilation (DA) system. This system was utilized to assimilate Himawari-8 aerosol optical thickness (AOT) into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to depict one SA air pollution outbreak event in spring 2018. Sensitivity tests for the assimilation system were conducted first to tune temporal localization lengths. Comparisons with independent Moderate Resolution Imaging Spectroradiometer (MODIS) and AErosol RObotic NETwork (AERONET) observations demonstrate that the AOT analysis and forecast fields have more reasonable diurnal variations by assimilating all the observations within a 12 h window, which are both better than assimilating the hourly observations in the current assimilation time slot. Assimilation of the entire window of observations with aerosol radiative effect activation significantly improves the prediction of downward solar radiation compared to the free-run experiment. The AOT assimilation with aerosol radiative effect activation led to a reduction in aerosol concentrations over SA, resulting in increased surface radiation, temperature, boundary layer height, and atmospheric instability. These changes facilitated air uplift, promoting aerosol transport from SA to the southeastern TP and leading to an increase in AOT in this region.

Including the Horizontal Observation Error Correlation in the Ensemble Kalman Filter: Idealized Experiments with NICAM-LETKF

Abstract Densely observed remote sensing data such as radars and satellites generally contain significant spatial error correlations. In data assimilation, the observation error covariance matrix is usually assumed to be diagonal, and the dense data are thinned or spatially averaged to compensate for neglecting the spatial observation error correlation. However, in theory, including the spatial observation error correlation in data assimilation can make better use of the dense data. This study performs perfect model observing system simulation experiments (OSSEs) using the nonhydrostatic icosahedral atmospheric model (NICAM) and the local ensemble transform Kalman filter (LETKF) to assess the impact of assimilating horizontally dense and error-correlated observations. The condition number of the observation error covariance matrix, defined as the ratio of the largest to smallest eigenvalues, is important for the numerical stability of the LETKF computation. A large condition number makes it difficult to compute the ensemble transform matrix correctly. Reducing the condition number by reconditioning is found effective for stable computation. The results show that including the horizontal observation error correlation with reconditioning makes the analysis more accurate but requires 6 times more computations than the case with the diagonal observation error covariance matrix. Significance Statement It is important to effectively utilize observations in data assimilation for more accurate weather prediction. Spatially dense observations are known to have an error correlation that is ignored in the data assimilation. This study explores assimilating dense observations by explicitly including observation error correlations with an idealized experiment. The results shows that the analysis is improved by including the observation error correlations. Also, the condition number of the observation error covariance matrix is essential for stable computations.

Improving the representation of the atmospheric boundary layer by direct assimilation of ground‐based microwave radiometer observations

AbstractIn a joint effort, MeteoSwiss and Deutscher Wetterdienst (DWD) address the need for improving the initial state of the atmospheric boundary layer (ABL) by exploiting ground‐based profiling observations that aim to fill the existing observational gap in the ABL. We implemented brightness‐temperature observations from ground‐based microwave radiometers (MWRs) in our data assimilation systems using a local ensemble transform Kalman filter (LETKF) with Radiative Transfer for TIROS Operational Vertical Sounder, ground‐based (RTTOV‐gb) as a forward operator. We were able to obtain a positive impact on the brightness temperature first guess and analysis, as well as a slight impact on the ABL humidity, using two MWRs at MeteoSwiss. These results led to a subsequent operational implementation of the observing system at MeteoSwiss. Furthermore, we performed an extensive set of assimilation experiments at DWD to investigate further various aspects such as the vertical localisation of selected single channels. We obtained a positive impact on the 6‐hr forecast of ABL temperature and humidity by assimilating two channels employing a dynamical localisation based on the sensitivity functions of RTTOV‐gb but also with a static localisation in a single‐channel setup. Our experiments indicate the importance of vertical localisation when using more than one channel, although reliable improvements are challenging to obtain without a larger number of observations for both assimilation and verification.

Continuous data assimilation of large eddy simulation by lattice Boltzmann method and local ensemble transform Kalman filter (LBM-LETKF)

We investigate the applicability of the data assimilation (DA) to large eddy simulations based on the lattice Boltzmann method (LBM). We carry out the observing system simulation experiment of a two-dimensional (2D) forced isotropic turbulence, and examine the DA accuracy of the nudging and the local ensemble transform Kalman filter (LETKF) with spatially sparse and noisy observation data of flow fields. The advantage of the LETKF is that it does not require computing spatial interpolation and/or an inverse problem between the macroscopic variables (the density and the pressure) and the velocity distribution function of the LBM, while the nudging introduces additional models for them. The numerical experiments with grids and 10% observation noise in the velocity showed that the root mean square error of the velocity in the LETKF with observation points ( of the total grids) and 64 ensemble members becomes smaller than the observation noise, while the nudging requires an order of magnitude larger number of observation points to achieve the same accuracy. Another advantage of the LETKF is that it well keeps the amplitude of the energy spectrum, while only the phase error becomes larger with more sparse observation. We also see that a lack of observation data in the LETKF produces a spurious energy injection in high wavenumber regimes, leading to numerical instability. Such numerical instability is known as the catastrophic filter divergence problem, which can be suppressed by increasing the number of ensemble members. From these results, it was shown that the LETKF enables robust and accurate DA for the 2D LBM with sparse and noisy observation data.

Comparative study of strongly and weakly coupled data assimilation with a global land–atmosphere coupled model

Abstract. This study explores coupled land–atmosphere data assimilation (DA) for improving weather and hydrological forecasts by assimilating soil moisture (SM) data. This study integrates a land DA component into a global atmospheric DA system of the Nonhydrostatic ICosahedral Atmospheric Model and the local ensemble transform Kalman filter (NICAM-LETKF) and performs both strongly and weakly coupled land–atmosphere DA experiments. We explore various types of coupled DA experiments by assimilating atmospheric observations and SM data simultaneously. The results show that analyzing atmospheric variables by assimilating SM data improves the SM analysis and forecasts and mitigates a warm bias in the lower troposphere where a dry SM bias exists. On the other hand, updating SM by assimilating atmospheric observations has detrimental impacts due to spurious error correlations between the atmospheric observations and land model variables. We also find that assimilating SM by strongly coupled DA is beneficial in the Sahel and equatorial Africa from May to October. These regions are characterized by seasonal variations in the precipitation patterns and benefit from updates in the atmospheric variables through SM DA during periods of increased precipitation. Additionally, these regions coincide with those identified in the previous studies, where a global initialization of SM would enhance the prediction skill of seasonal precipitation.

Assessment of Inflation Schemes on Parameter Estimation and Their Application in ENSO Prediction in an OSSE Framework

The ensemble Kalman filter is often used in parameter estimation, which plays an essential role in reducing model errors. However, filter divergence is often encountered in an estimation process, resulting in the convergence of parameters to the improper value and finally in parameter estimation failure. To alleviate this degeneration, various covariance inflation schemes have been proposed. In this study, I examined six currently used inflation schemes: fixed inflation, conditional covariance inflation, modified estimated parameter ensemble spread, relaxation-to-prior perturbations, relaxation-to-prior spread, and new conditional covariance inflation. The six schemes were thoroughly explored using the Zebiak–Cane model and the local ensemble transform Kalman filter in the observing system simulation experiment framework. Emphasis was placed on the comparison of these schemes when it came to estimating single and multiple parameters in terms of oceanic analyses and resultant El Niño–Southern Oscillation (ENSO) predictions. The results showed that the new conditional covariance inflation scheme had the best results in terms of the estimated parameters, resultant state analyses, and ENSO predictions. In addition, the results suggested that better parameter estimation yields better state simulations, resulting in improved predictions. Overall, this study provides viable information for selecting inflation schemes for parameter estimation, offering theoretical guidance for constructing operational assimilation systems.

Impact of Gaussian Transformation on Cloud Cover Data Assimilation for Historical Weather Reconstruction

Abstract Old descriptive diaries are important sources of daily weather conditions before modern instrumental measurements were available. A previous study demonstrated the potential of reconstructing historical weather at a high temporal resolution by assimilating cloud cover converted from descriptive diaries. However, cloud cover often exhibits a non-Gaussian distribution, which violates the basic assumptions of most data assimilation schemes. In this study, we applied a Gaussian transformation (GT) approach to cloud cover data assimilation and conducted observing system simulation experiments (OSSEs) using 20 observation points over Japan. We performed experiments to assimilate cloud cover with large observational errors using the Global Spectral Model (GSM) and a local ensemble transform Kalman filter (LETKF). Without GT, meridional wind and temperature exhibited deteriorations in the lower troposphere compared with the experiment with no observations. In contrast, GT reduced the 2-month root-mean-square errors (RMSEs) by 5%–15% throughout the troposphere for wind, temperature, and specific humidity fields. Significant improvements include zonal wind at 500 hPa and temperature at 850 hPa with 6.4% and 7.3% improvements by GT, respectively, compared with the experiment without GT. We further demonstrate that the additional GT application to the precipitation background field improves precipitation estimation by 12.2%, with pronounced improvements over regions with monthly precipitation of less than 150 mm. We also explored the impact of cloud cover GT on a global scale and confirmed improvements extending from around the observation sites. Our results demonstrate the potential of GT in high-resolution historical weather reconstruction using old descriptive diaries. Significance Statement To reconstruct the historical weather, cloud cover information from old diaries can be used by incorporating high-resolution model simulations. However, cloud cover is not normally distributed and violates an important assumption when combining cloud cover observations with model simulations. Our results demonstrate that transforming the cloud cover distribution into a normal distribution could improve wind speed, temperature, and humidity fields in the model. We demonstrate the critical role of the transformation in a nonnormally distributed variable when combined with models and show the potential of diary-based weather information to reconstruct historical weather.

Assessing the Ensemble Assimilation System for Forecasting the Genesis of Tropical Cyclone Sonca (2017) in the Vietnam East Sea using a Lagrangian Approach

: With the aim of evaluating the cyclogenesis predictability of Sonca (2017) using an ensemble prediction system from the local ensemble transform Kalman filter (LETKF) to model WRF, the study proposed a dynamical assessment method with a Lagrangian approach that is stable and independent of time and reference systes. With this approach, the study can clearly distinguish between cases of tropical cyclone (TC) formation and non-formation in the ensemble forecast members. By evaluating the thermodynamic conditions of the low-level circulation between the forming and non-forming groups, the study found that the ocean-atmosphere interaction and the vortex merger process play imnportant roles in the formation of Sonca. The evaluation of the formation location showed that the detected vortex centers from the forming members tended to concentrate near the monsoon trough axis, indicating that this is a highly favorable region for vortex formation in the Vietnam East Sea.
  
  
  

Reconstructing cardiac electrical excitations from optical mapping recordings.

The reconstruction of electrical excitation patterns through the unobserved depth of the tissue is essential to realizing the potential of computational models in cardiac medicine. We have utilized experimental optical-mapping recordings of cardiac electrical excitation on the epicardial and endocardial surfaces of a canine ventricle as observations directing a local ensemble transform Kalman filter data assimilation scheme. We demonstrate that the inclusion of explicit information about the stimulation protocol can marginally improve the confidence of the ensemble reconstruction and the reliability of the assimilation over time. Likewise, we consider the efficacy of stochastic modeling additions to the assimilation scheme in the context of experimentally derived observation sets. Approximation error is addressed at both the observation and modeling stages through the uncertainty of observations and the specification of the model used in the assimilation ensemble. We find that perturbative modifications to the observations have marginal to deleterious effects on the accuracy and robustness of the state reconstruction. Furthermore, we find that incorporating additional information from the observations into the model itself (in the case of stimulus and stochastic currents) has a marginal improvement on the reconstruction accuracy over a fully autonomous model, while complicating the model itself and thus introducing potential for new types of model errors. That the inclusion of explicit modeling information has negligible to negative effects on the reconstruction implies the need for new avenues for optimization of data assimilation schemes applied to cardiac electrical excitation.

Improvements of Lagrangian Data Assimilation Tested in the Gulf of Mexico

Abstract This study extends initial work by Sun and Penny and Sun et al. to explore the inclusion of path information from surface drifters using an augmented-state Lagrangian data assimilation based on the local ensemble transform Kalman filter (LETKF-LaDA) with vertical localization to improve analysis of the ocean. The region of interest is the Gulf of Mexico during the passage of Hurricane Isaac in the summer of 2012. Results from experiments with a regional ocean model at eddy-permitting and eddy-resolving model resolutions are used to quantify improvements to the analysis of sea surface velocity, sea surface temperature, and sea surface height in a data assimilation system. The data assimilation system assimilates surface drifter positions, as well as vertical profiles of temperature and salinity. Data were used from drifters deployed as a part of the Grand Lagrangian Deployment beginning 20 July 2012. Comparison of experiment results shows that at both eddy-permitting and eddy-resolving horizontal resolutions Lagrangian assimilation of drifter positions significantly improves analysis of the ocean state responding to hurricane conditions. These results, which should be applicable to other tropical oceans such as the Bay of Bengal, open new avenues for estimating ocean initial conditions to improve tropical cyclone forecasting.

Assimilating GNSS TEC with an LETKF over Yunnan, China

A robust ionospheric model is indispensable for providing the atmospheric delay corrections for global navigation satellite system (GNSS) navigation and positioning and forecasting the space environment. The accuracy of ionospheric models is limited due to the simplified model structures. Complicated spatiotemporal variations in total electron content (TEC) biases between GNSS and international reference ionosphere (IRI) suggest a robust strategy to optimally combine GNSS and IRI TEC for high-precision modeling. In this paper, we propose a novel ionospheric data assimilation method, which is a local ensemble transform Kalman filter (LETKF), to construct an ionospheric model over Yunnan in southwestern China. We used the LETKF method to assimilate the ionospheric TEC extracted from GNSS observations in Yunnan into the IRI-2016 model. The experimental results indicate that the ionospheric data assimilation has a more pronounced improvement effect on the IRI empirical model during periods of geomagnetic quiet than during periods of geomagnetic disturbance. On quiet magnetic days, the skill score (SKS) of the assimilation is 0.60 and the root mean square error (RMSE) values before and after assimilation are 5.08 TECU and 2.02 TECU, respectively. The correlation coefficient after assimilation increases from 0.94 to 0.99. On magnetic storm days, the SKS of the assimilation is 0.42 and the RMSE values before and after assimilation are 5.99 TECU and 3.46 TECU, respectively. The correlation coefficient after assimilation increases from 0.98 to 0.99. The results suggest that the LETKF algorithm can be considered an effective method for ionospheric data assimilation.

The Local Unscented Transform Kalman Filter for the Weather Research and Forecasting Model

In this study, the local unscented transform Kalman filter (LUTKF) proposed in the previous study estimates the state of the Weather Research and Forecasting (WRF) model through local analysis. Real observations are assimilated to investigate the analysis performance of the WRF-LUTKF system. The WRF model as a regional numerical weather prediction (NWP) model is widely used to explain the atmospheric state for mesoscale meteorological fields, such as operational forecasting and atmospheric research applications. For the LUTKF based on the sigma-point Kalman filter (SPKF), the state of the nonlinear system is estimated by propagating ensemble members through the unscented transformation (UT) without making any linearization assumptions for nonlinear models. The main objective of this study is to examine the feasibility of mesoscale data assimilations for the LUTKF algorithm using the WRF model and real observations. Similar to the local ensemble transform Kalman filter (LETKF), by suppressing the impact of distant observations on model state variables through localization schemes, the LUTKF can eliminate spurious long-distance correlations in the background covariance, which are induced by the sampling error due to the finite ensemble size; therefore, the LUTKF used in the WRF-LUTKF system can efficiently execute the data assimilation with a small ensemble size. Data assimilation test results demonstrate that the LUTKF can provide reliable analysis performance in estimating the WRF model state with real observations. Experiments with various ensemble size show that the LETKF can provide better estimation results with a larger ensemble size, while the LUTKF can achieve accurate and reliable assimilation results even with a smaller ensemble size.

Implementation and Testing of Radar Data Assimilation Capabilities Within the Joint Effort for Data Assimilation Integration Framework With Ensemble Transformation Kalman Filter Coupled With FV3‐LAM Model

AbstractCapabilities to directly assimilate radar data are implemented within the local ensemble transform Kalman filter (LETKF) and the gain‐form LETKF (LGETKF) algorithms of the Joint Effort for Data assimilation Integration (JEDI) system. The capabilities are evaluated for the analysis and forecast of a severe convection case of 20 May 2019 in the Southern Great Plains using the limited area model version of the FV3 dynamical core (FV3‐LAM) from a recent release for Short‐Range Weather Application (SRW App). The LETKF and LGETKF implementations are shown to produce analyses and short‐range forecasts comparable to those using the ensemble square‐root Kalman Filter (EnSRF) within the Gridpoint Statistical Interpolation (GSI) framework used by current NCEP operational models. In addition, LGETKF retaining only 60% variances for model‐space vertical localization performs similarly to LGETKF retaining 99% of variance and LETKF using observation error‐based vertical localization. JEDI LETKF shows better parallel scalability than LGETKF and GSI EnSRF.

Estimation of AMSU-A Radiance Observation Impacts in an LETKF-Based Atmospheric Global Data Assimilation System: Comparison with EFSO and Observing System Experiments

Abstract This work assesses the contribution of assimilating AMSU-A satellite-based radiance measurements to a global data assimilation system based on an atmospheric general circulation model and the local ensemble transform Kalman filter (LETKF). The radiance measurements were from three channels that are sensitive to the upper troposphere and lower stratosphere. The contribution of these measurements, or AMSU-A observation impact, was estimated both through ensemble-based forecast sensitivity to observations (EFSO) and observing system experiments (OSEs). Two streams of data-denial experiments for the AMSU-A observations were performed for about one month during winter in each hemisphere. The OSEs quantified the accumulated observation impact by cycling (repeating) data denials: including AMSU-A observations reduced the total observation impact for all observations of each data assimilation cycle. In contrast, EFSO estimated AMSU-A to increase the total observation impact. The opposing effects were attributed to the accumulated observation impact in the OSEs; the accumulation could stabilize the data assimilation cycles. In both experiments, the accumulated observation impact of AMSU-A was strongest in the upper troposphere, particularly in the austral midlatitudes where westerly jets exist and observations of other types are sparse. EFSO also assessed AMSU-A to have the most beneficial observation impact in similar locations. The AMSU-A observation impact tended to accumulate just downstream of where EFSO estimated the beneficial observation impact signals. The accumulated AMSU-A observation impact was tied to dynamic processes in the upper-tropospheric and general stratospheric circulation. Therefore, EFSO helps estimate the beneficial distributions of AMSU-A accumulated observation impact by considering their dynamical propagation. Significance Statement The Advanced Microwave Sounding Unit-A (AMSU-A) satellite radiance assimilation technique was successfully integrated into the Atmospheric General Circulation Model for the Earth Simulator (AFES)–LETKF data assimilation system. We conducted OSEs and used EFSO to assess the AMSU-A observation impact. The two estimation methods identified opposite observation impacts due to the cycling (repeating) OSEs of the AMSU-A observations. We interpreted the causes of the opposite estimations. However, even for the cycling OSEs, EFSO appeared to help estimate distributions of the accumulated observation impact. It is important to consider the dynamical propagation of accumulated observation impact in general circulation.

Ensemble-based data assimilation of significant wave height from Sofar Spotters and satellite altimeters with a global operational wave model

An ensemble-based method for wave data assimilation is implemented using significant wave height observations from the globally distributed network of Sofar Spotter buoys and satellite altimeters. The Local Ensemble Transform Kalman Filter (LETKF) method generates skillful analysis fields resulting in reduced forecast errors out to 2.5 days when used as initial conditions in a cycled wave data assimilation system. The LETKF method provides more physically realistic model state updates that better reflect the underlying sea state dynamics and uncertainty compared to methods such as optimal interpolation. Skill assessment far from any included observations and inspection of specific storm events highlight the advantages of LETKF over an optimal interpolation method for data assimilation. This advancement has immediate value in improving predictions of the sea state and, more broadly, enabling future coupled data assimilation and utilization of global surface observations across domains (atmosphere-wave-ocean).